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Abstract 

We present a new limit on the production of a light dark-force mediator with 
the KLOE detector at DA<£NE. This boson, called U, has been searched for in the 
decay (f> — > rjU, with the hypothesis U — > e + e~ , analyzing the decay n — > 7r°7r 7r° 
in a data sample of 1.7 fb _1 integrated luminosity. This search is combined with a 
previous result obtained using the decay n — > tt + tt~ tt° , increasing substantially the 
sensitivity. No structures are observed in the e + e~ invariant mass distribution over 
the background. Combining the two n decay modes, we set an upper limit at 90% 
C.L. on the ratio between the U boson coupling constant and the fine structure 
constant of a' /a < 1.5 x 1(T 5 for 30 < M v < 420 MeV and a' /a < 5.0 x 10~ 6 
for the sub-region 60 < Mjj < 190 MeV. This result assumes the Vector Meson 
Dominance expectations for the (pnj* transition form factor. The dependence of 
this limit on the transition form factor has also been studied. 

Key words: e + e~ collisions, dark forces, gauge vector boson 
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1 Introduction 



There are theories beyond the Standard Model (SM) that postulate the ex- 
istence of new bosons, mediators of some hidden gauge group, under which 
ordinary matter is uncharged [Tf2] . Such theories have been recently invoked 
to account for some puzzling astrophysical observations and predict under the 
hypothesis that at least one of the new bosons is neutral, relatively light and 
weakly coupled with ordinary matter [3~P] . The coupling is represented by 
adding to the Lagrangian a kinetic mixing term between the new (U) boson 
and the photon, whose strength is measured by a parameter e that could be as 
large as ~ 10 -3 . Due to this mixing, the U boson can be produced in processes 
involving SM particles, and can decay into e + e~, 7r + 7r~..., depending 

on its mass and on model-specific details. 

Searches for the U boson have been recently performed at e-p fixed target 
facilities [HE], and, in conjunction with a dark Higgs, at e + e~ colliders [8], 
with null results. The U can also be searched for in vector (V) to pseudoscalar 
(P) meson decays, with a rate that is e 2 times suppressed with respect to 
the ordinary V —> P7 transitions [9]. Since the U is supposed to decay to 
e + e~ with a non negligible branching ratio, V — > PU events will produce 
a sharp peak in the invariant mass distribution of the electron-positron pair 
over the continuum background due to Dalitz decay events V — > Pe + e~ . Using 
this approach, KLOE has already published a limit on the existence of the U 
boson [10], studying <fi — > 7]e + e~ decays, where the 77 meson was tagged by 
its 7r + 7T _ 7r° decay. In this letter, we present an update of this analysis, which 
improves background rejection for the already used data sample (1.5 fb _1 ) 
and increases statistics by a factor about three, exploiting also the neutral 
f] — > 7r 7r 7T decay chain, with a data sample of 1.7 fb _1 . 



2 The KLOE detector 



DA$NE , the Frascati 0-factory, is an e + e~ collider running at center-of-mass 
energy of ~ 1020 MeV. Positron and electron beams collide at an angle of 7T-25 
mrad, producing <ft mesons nearly at rest. The KLOE experiment operated at 
this collider from 2000 to 2006, collecting 2.5 fb -1 . The KLOE detector consists 
of a large cylindrical Drift Chamber (DC), surrounded by a lead-scintillating 
fiber electromagnetic calorimeter (EMC), all embedded inside a superconduct- 
ing coil, providing a 0.52 T axial field. The beam pipe at the interaction region 
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is a sphere with 10 cm radius, made of a 0.5 mm thick Beryllium- Aluminum 
alloy. The drift chamber 4 m in diameter and 3.3 m long, has 12,582 all- 
stereo tungsten sense wires and 37,746 aluminum field wires, with a shell made 
of carbon fiber-epoxy composite with an internal wall of ~ 1 mm thickness. 
The gas used is a 90% helium, 10% isobutane mixture. The momentum resolu- 
tion is a(p±)/p± ~ 0.4%. Vertices are reconstructed with a spatial resolution of 
~ 3 mm. The calorimeter [T2] , with a readout granularity of ~ (4.4 x 4.4) cm 2 , 
for a total of 2440 cells arranged in five layers, covers 98% of the solid angle. 
Each cell is read out at both ends by photomultipliers, both in amplitude and 
time. The energy deposits are obtained from the signal amplitude while the 
arrival times and the particles positions are obtained from the time differences. 
Cells close in time and space are grouped into energy clusters. Energy and time 
resolutions are a E /E = 5.7%/^E (GeV) and a t = 57 ps/^E (GeV) © 100 ps, 
respectively. The trigger [13] uses both calorimeter and chamber information. 
In this analysis the events are selected by the calorimeter trigger, requiring 
two energy deposits with E > 50 MeV for the barrel and E > 150 MeV for the 
endcaps. Data are then analyzed by an event classification filter [H], which 
selects and streams various categories of events in different output files. 



3 Event selection 

To improve the search for the U boson, we have carried out the analysis of 
the process <f> —> r]U, U — > e + e~, adding the decay channel 77 — » 7r 7r°7r°, to 
the previously used, rj — > 7r + 7r~7r°. The new search has been performed on a 
data sample of 1.7 fb _1 , corresponding approximately to 6 x 10 9 produced <f> 
mesons. The MC simulation for the — > rjU decay has been developed ac- 
cording to [9], with a flat distribution in the e + e~ invariant mass, M ee , while 
the irreducible background <fi — > r/e + e~, r] — > mm, has been simulated accord- 
ing to a Vector Meson Dominance parametrization [15] . All MC productions, 
including all other cf> decays, take into account changes in DA$NE operation 
and background conditions on a run-by-run basis. Corrections for data-MC 
discrepancies in cluster energies and tracking efficiency, evaluated with radia- 
tive Bhabha scattering and — ?• p7r event samples, respectively, have been 
applied. 

As first analysis step for the neutral r\ decay channel, a preselection is per- 
formed requiring: 

(1) two opposite charge tracks with point of closest approach to the beam 
line inside a cylinder around the interaction point (IP), of 4 cm transverse 
radius and 20 cm length; 

(2) six prompt photon candidates i.e. energy clusters with E > 7 MeV not 
associated to any track, in an angular acceptance | cos 6^1 < 0.92 and in 
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Fig. 1. Recoiling mass against the e + e~ pair for the data sample after preselection 
cuts. The <f> — y Tj e + e~ signal is clearly visible as the peak corresponding to the rj 
mass. 

the expected time window for a photon (|T 7 — Ry/c\ < MIN(3<r i , 2 ns)); 
(3) a loose cut on the six photon invariant mass: 400 < M 67 < 700 MeV. 

After this selection, a peak corresponding to the r\ mass is clearly observed in 
the distribution of the recoil mass against the e + e~ pair, M recoi i(ee) (Fig. [Q. 
The second peak at ~ 590 MeV is due to Kg —> 7i + tt~ decays with wrong 
mass assignment. To select — > r]e + e~ events, a 3a cut is applied on this 
variable, 536.5 < M recoil (ee) < 554.5 MeV. The retained sample has ~ 20% 
residual background contamination, constituted by <fi — > 777, —> K S K L and 
e + e" — > ojir (about 50%, 35% and 15% of the whole background contribution, 
respectively). In Fig. [2], the comparison between data and Monte Carlo (MC) 
events for the M ee and cos ^* distributions is shown at this analysis level. The 
^F* variable is the angle between the directions of the rj and the e + in the e + e~ 
rest frame. Photon conversion events are concentrated at M ee ~ 30 MeV and 
cos^* < 0.6, while the other backgrounds cover the M ee > 400 MeV region 
and are uniformly distributed in cos\l/*. 

The — > 777 background contamination is mainly due to events where a photon 
converts to an e + e~ pair on the beam pipe (BP) or drift chamber walls (DCW). 
After tracing back the tracks of the two e + /e~ candidates, these events are 
efficiently rejected by reconstructing the invariant mass (M ee ) and the distance 
(D ee ) of the track pair both at the BP and DCW surfaces. Both variables are 
expected to be small for photon conversion events, so that this background is 
removed by rejecting events with: [ M ee (BP) < 10 MeV and D ee (BP) < 2 
cm ] or [ M ee (DCW) < 120 MeV and D ee (DCW) < 4 cm ]. 

At this stage of the analysis, the surviving background is dominated by events 
with two charged pions in the final state, and it is rejected by exploiting the 
timing capabilities of the calorimeter. When an energy cluster is associated 
to a track, the time-of-flight (ToF) to the calorimeter is evaluated both using 
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Fig. 2. (j) — >• r\ e + e , rj — > 7r 7r°7r° events: data-MC comparison for M ee (top) and 
cos^* distributions (bottom) after the M recoi i(ee) cut. 



the track trajectory (T track = L track / (3c) and the calorimeter timing (T cluster ). 
The AT = T track — T c i ustcr variable is then evaluated in the electron hypothesis 
(AT e ). In order to be fully efficient on signal, events with either an e + or an 
e~ candidate inside a 3cr window around AT e = are retained for further 
analysis. 

At the end of the analysis chain, 30577 events are selected, with ~ 3% back- 
ground contamination (Fig. [3]). The analysis efficiency, defined as the ratio 
between the number of events surviving analysis cuts and that of all gener- 
ated events, is ~ 15% for low e + e~ invariant masses, increasing up to 30% for 
higher M ee values. 

The analysis of the decay channel r] — > 7r + 7r~7r° is the same as described in 
[TU] . with the addition of a cut on the recoil mass to the e + e~7r + 7r~ system, 
which is expected to be equal to the ir° mass for signal events. In Fig. H] top, 
data-MC comparison shows some residual background contamination in the 
tails of the distribution, which are not well described by our simulation. A cut 
100 < M recoi i(ee7r7r) < 160 MeV is then applied. The effect of this cut on the 
M ee variable is shown in Fig. H] bottom. 
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Fig. 3. 4> — > r]e + e~, rj — > 7r 7r°7r° events: data-MC comparison for M ee (top) and 
cos ^* distributions (bottom) at the end of the analysis chain. 

4 Upper limit evaluation on U boson production 



The evaluation of the upper limit on the U boson production in the <fi ~^ vU 
process is obtained combining the two rj decay channels. The resolution of 
the e + e~ invariant mass is less than 2 MeV over the whole M ee range. The 
determination of the limit is done by varying the My mass, with 1 MeV 
step, in the range between 5 and 470 MeV. Only five bins (5 MeV width) 
of the reconstructed M ee variable, centered at Mjj are considered. For each 
channel, the irreducible background, b(Mjj), is extracted directly from our data 
after applying a bin-by-bin subtraction of the non-irreducible backgrounds 
and correcting for the analysis efficiency. The M ee distribution is then fit, 
excluding the bins used for the upper limit evaluation. The parametrization 
of the fitting function has been taken from Ref. [15]. The (prj^f* transition form 
factor is parametrized as 

F M) = I _ l q2/K2 , (1) 

with q = M ee . Free parameters are A and a normalization factor. The spread 
of the extracted parameters is contained within the statistical error of the fit 
done on the whole M ee mass range, shown in Fig. El as expected from the 
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Fig. 4. cj) — > rje + e~, t] — > tt + tt^tt° analysis. Top: data-MC comparison for the recoil 
mass against the e + e~ir + 7r~ system. Bottom: M ee distribution before (open circles) 
and after (black dots) the cut on M reco ii(ee7T7r) 
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Fig. 5. Fit to the corrected M ee spectrum for the Dalitz decays 4> —> r]e + e~, with 
7] — > 7T 7r 7r° (top) and r] — > ir + ir~Tr° (bottom). 

overall good description of the M ee shape for both rj decay channels. 

The exclusion limit on the number of events for the <f> — > r\ U signal as a 
function of My is obtained with the CLs technique [16] , with the same method 
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Fig. 6. Top: upper limit at 90% C.L. on number of events for the decay chain 
4> — > T) U, U — > e + e~ , with r\ — > 7r 7r 7r° and r\ — > 7r + 7r~7r . Bottom: smoothed upper 
limit at 90% C.L. on BR(c/> — > r/U — > r/e + e~), obtained separately for the two r/ 
decay channels and from the combined analysis. 

used in our previous analysis [10] . The limit is extracted both for each t] decay 
channel and in a combined way. The combination uses simultaneously the two 
samples weighted with luminosity, efficiency and relative branching ratios. 
The systematic error on the background knowledge Ab(M ee ) is evaluated, for 
each M[/ value, changing by one standard deviation the two fit parameters. A 
Gaussian spread of width A6(M ee ) on the background distribution is applied 
while evaluating CLs. In Fig. |6] top, the upper limit at 90% C.L. on the 
number of events for the decay chain — > 77 U, U — > e + e~, is shown for both 
f] — > 7r 7r°7r° and rj — > 7r + 7r~7r°, separately evaluated. In Fig. [6] bottom, the 
smoothed upper limit on the branching fraction for the process <f> — > r\ U, 
U — > e + e~, obtained from the combined method is compared with evaluations 
from each of the two decay channels. 

The exclusion plot in the a' /a = e 2 vs Mjj plane, where a' is the coupling of 
the U boson to electrons and a is the fine structure constant, has been finally 
derived assuming the relation [9]: 

a(e+e~ -+ ^ r,U) = e 2 |i^«)| 2 a(e+ e - 0^ m) , 

| (2) 

with X(m 2 , m 2 , m^) = [1 + m 2 / (m\ — m\)] 2 — Amfm 2 / (mf — m^) 2 . We assume 
that the U boson decays only to lepton pairs, with equal coupling to e + e~ and 
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The extraction of the limit on the a'/ a parameter is related to the parametriza- 
tion of the form factor (Eq. [2]), and thus to the A parameter in Eq. [TJ The 
SND experiment measured the form factor slope, = dF/dq 2 \ q2=0 = A~ 2 , 
obtaining b^ v = (3.8 ± 1.8) GeV~ 2 , whith a central value different from theo- 
retical predictions based on VMD (b^ ~ 1 GeV -2 ) [5], although in agreement 
within the error. In Fig. [7] the smoothed exclusion plot at 90% CL. on a'/ a 
is compared with existing limits in the same region of interest [6. 71TTT] . The 
evaluation is done using both the experimental and the theoretical values of 
the form factor slope. The two resulting curves overlap at low M ee values, 
while the limit obtained using the SND measurement gives an increasingly 
larger exclusion region up to ~ 400 MeV, rejoining the other curve at the end 
of the phase space. Having the experimental value of b^ an uncertainty of 
~ 50%, we conservatively use the curve obtained with theoretical predictions, 
resulting in a limit of: a'/a < 1.5 x 10~ 5 for 30 < Mu < 420 MeV, and even 
better for the sub-region 60 < My < 190 MeV: a'/a < 5.0 x 10 -6 . Comparing 
our result with the previous KLOE measurement, reported as dotted line in 
Fig. we improve the upper limit of about a factor two, both including the 
7] — > 7r°7r°7r final state and improving the analysis for the i] — > 7r + 7r~7r° decay 
channel. This result reduces the region of the U boson parameters that could 
explain the observed discrepancy between measurement and Standard Model 
prediction of the muon anomalous magnetic moment, a^, ruling out masses in 
the range 50-460 MeV. 
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